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Abstract 



The cross sections of heavy quark production in pp collision and 
for their photo- and electroproduction are calculated in the framework 
CN | of QCD. The virtual nature of the interacting gluons as well as their 

transverse motion and different polarizations are taken into account. 



£f~) ■ The obtained cross sections exhibit more rapid growth with the initial 

energy than the parton model predictions and the pt distributions of 
■ produced heavy quarks are more smooth. 
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1 INTRODUCTION 



The investigation of the production of heavy quarks in high energy hadron 
processes provides a method for studying the internal structure of hadrons. 
Realistic estimates of the cross section of the heavy quark production are 
necessary in order to plan experiments on existing and future accelerators. 
The predicted cross section values are usually obtained in the parton model 
framework and depend significantly on the quark and gluon structure func- 
tions. These functions are practically unknown experimentally for very small 
values of Bjorken variable x. However it is just the region that dominates in 
the heavy quark production at high energies since the main contribution to 
heavy quark hadroproduction cross section comes from the x values of the 
order of x ~ 2m T /^/s , where ttit = \J m 2 + -p\ is transverse mass of the 
produced quark. 

The Gribov-Lipatov-Altarelli-Parisi (GLAP) evolution equation is used 
usually to calculate the structure functions. It sums up in leading loga- 
rithm approximation (LLA) all the QCD diagram contributions proportional 
to [a s In q 2 ) n but it does not take into account the terms proportional to 
(a s ln l/x) n , so this approximation does not give the correct asymptotic be- 
haviour of the structure function in the small x region. For the correct 
description of these phenomena not only the terms of the form (a s In q 2 ) n 
have to be collected in the Feynman diagrams but also the terms (a s In l/x) n 
and (a s ln q 2 ■ In l/x) n . 

Another problem that appears at x ~ is that of the absorption (screen- 
ing) corrections which must stop the increase of the cross section when x — > 
in accordance with the unitarity condition. It can be interpreted as satura- 
tion of the parton density. For relatively small virtuality q 2 < q 2 (x) the gluon 
structure function behaves as xG(x, q 2 ) ~ q 2 R 2 , so that the cross section for 
the interaction of a point-like parton with the target, a ~ (l/q 2 )xG(x, q 2 ) ~ 
R 2 , obeys the unitarity condition. The quantity q 2 (x) can be treated as a 
new infrared-cutoff parameter which plays the role of a typical transverse 
momentum of partons in the parton cascade of the hadron in semihard pro- 
cesses. The behaviour of qo(x) was discussed in ref. [6]. It increases with 
log(l/:r) and at x = 0.01 — 0.001 the values of qo(x) are about 2-4 GeV. 

The predictions for the cross section value of heavy quark pair produc- 
tion are based usually on the parton model calculations [1, 2]. In this model 
all particles involved are assumed to be on mass shell, having only longitu- 
dinal component of the momentum and the cross section is averaged over 
two transverse polarizations of the gluons. The virtualities q 2 of the ini- 
tial partons are taken into account only through their densities. The latter 
are calculated in LLA through GLAB evolution equation. The probabilis- 
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tic picture of noninteracting partons underlies this way of proceeding. In 
the region where the transverse mass tut of the produced heavy quark Q is 
close to qo(x) the dependence of the amplitude of the subprocesses gg — ► QQ 
or 7g — > QQ on the virtualities and polarizations of the gluons should be 
taken into account, i.e., the matrix elements of these subprocesses should be 
calculated more accurately than is usually done in the parton model. The 
matrix elements of the QCD subprocesses accounting for the virtualities and 
polarizations of the gluons are very complicated. We present here only the 
results for main and simplest subprocess in each case, gg — ► QQ (~ a 2 ) for 
hadroproduction and jg — ► QQ (~ a s ) for photo- and electroproduction. 
The contributions of high-order subprocesses can be essential but our aim 
is to discuss the qualitative difference between our results and the parton 
model predictions that can be done on the level of low-order diagrams. 

In this paper we calculate the cross section for three kinds of processes: for 
the photoproduction of heavy quarks, for electroproduction (in deep inelastic 
scattering) and for their production in hadron-hadron collisions. The latter 
case was concidered in our previous paper || , however the structure function 
we used there was too rough numerically although it reproduced properly all 
the qualitative features of the semihard interactions. In the present paper we 
use two gluon structure functions. The first one is Morfm-Tung (MT S-DIS) 
M parametrization which fits experimental data with high accuracy. The 
second was obtained in || using the same MT S-DIS boundary condition 
(at Q 2 = Qq) with accounting for the shadow corrections. We correct also a 
numerically small error in the matrix element presented in ||. 



2 CROSS SECTIONS OF HEAVY FLAVOUR 
PRODUCTION IN QCD 



The cross section of heavy quarks photo and electroproduction is given 
schematically by the graphs in Fig.l while Fig.2 presents the sum of dia- 
grams for their production in pp collision. The main contribution to the 
cross section at small x is known to come from gluons. The lower and upper 
(for pp case) ladder blocks present the two-dimensional gluon distribution 
ip(x, qf) and <f(x, q 2 ) which are the functions of the fraction (x and y) of the 
longitudinal momentum of the initial hadron and the gluon virtuality. Their 
distribution over x and transverse momenta qx in hadron is given in semihard 
theory by function (p(x, q 2 ). It differs from the usual function G(x, q 2 ): 




(1) 
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Such definition of <p(x, q 2 ) makes possible to treat correctly the effects arising 
from gluons virtualities. The exact expression for this function can be ob- 
tained as a solution of the evolution equation which, contrary to the parton 
model case, is nonlinear due to interactions between the partons in small x 
region. 

In what follows we use Sudakov decomposition for quarks' momenta pi j2 
through the momenta of colliding hadrons pa and ps {jp\ — P% — 0) an d 
transverse ones Pi j2 t- 

Pl,2 = x 1>2 Pb + UijPa + Pl,2T- (2) 

The differential cross sections of heavy quarks photo-, electro-Q and hadropro- 
duction have the form:| 



dcrph oi em eQ 



dyldy%d 2 p 1T d 2 p 2T (2vr) 4 (s 



4^ J d 2 q 2T 5(q 2T ~ Pit - P2T)S(yi + y 2 -l) 



x ^v(qlx)\M ph \ 2 ± (3) 

<?2 Vl 



d<J e i,T,L o> em eQ 



dytdy* 2 d 2 p 1T d 2 p 2T (2vr) 4 (s 



7p J d 2 q 2T 5(q 2T - Pit - P2t) <%i +yz-l) 



x ^^qlx)\M el>T>L \ 2 -, (1) 



12 V2 



d<Jpp 



1 1 f 

~]8~T)2 J d2 ^Td 2 q 2 T5(qiT + q2T - Pit - P2t) 



dy\dy* 2 d 2 p lT d 2 p 2T (2tt 

x as{ f )as{ f\ {qly)^qlx)\M pp \ 2 . (5) 

Here s = IpaPb , <Zi,2T are the gluons' transverse momenta and y\ 2 are the 
quarks' rapidities in the photon-, 7*- or hadron-hadron c.m.s. frames, 

Xi= m^ e - y ^ X2= m^ e - y ^ x = Xi + X2 

yi =^ e yi, y 2 = m^ e y^ y = y 1 + y 2 . 

\Mphl 2 , \M e i t T,L,\ 2 and |M pp | 2 are the squares of the matrix elements for 
the cases of photo-, electro- and hadroproduction, respectively, and indices 
T and L denotes the transverse and longitudinal parts. 

1 &ei,T,L denotes the virtual photon-proton transverse and longitudinal cross sections. 
2 We put the argument of as to be equal to gluon virtuality, which is very close to the 
BLM scheme @; (see also §]). 
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In LLA kinematic 



qi - VPa + qiT, <?2 - xp B + qiT- (7) 

so 

q\ - -5iT» 4 - -<&t- (8) 

2 2 

(The more accurate relations are q\ = — fq^ q\ — but we are working 

in the kinematics where x, y ~ 0). 

The matrix elements M are calculated in the Born order of QCD without 
standart simplifications of the parton model since in the small x domain there 
are no grounds for neglecting the transverse momenta of the gluons qyr and 
q2T in comparision with the quark mass and the parameter qo(x). In the axial 
gauge p B A^ = the gluon propagator takes the form D^q) = d fll/ (q)/q 2 , 

d, v {q) = 8, v - {ftft + q u p%)/(pBq). (9) 

For the gluons in t— channel the main contribution comes from the so called 
'nonsense' polarization , which can be picked out by decomposing the 
numerator into longitudinal and transverse parts: 

M?) = 2 (PbPa + PaPb)/s + Si « 2p BP y S = g; u . (10) 

The other contributions are suppressed by the powers of s. Since the sums of 
the diagrams in Fig.l and Fig. 2 are gauge invariant in the LLA, the transver- 
sality condition for the ends of gluon line enables one to replace p\ by —qi T /x 
in expression for gr" . Thus we get 

d, v {q) « -2^ (11) 
x s 

or 

d, v {q) « 1— (12) 
xys 

if we do such a trick for vector ps too. Both these equations for d^ v can be 
used but for the form (|11|) one has to modify the gluon vertex slightly (to 
account for several ways of gluon emission — see ref . ) : 

2 

Kff = — l( x y s - qIt) Qit - q 2 \Tq v 2T + 2x (qirq2T) p v b \- ( 13 ) 

As a result the colliding gluons can be treated as aligned ones and their po- 
larization vectors are directed along the transverse momenta. Ultimately, the 
nontrivial azimuthal correlations must arise between the transverse momenta 
Pit and p 2 T of the heavy quarks. 
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In the case of photoproduction the alignment of the gluon makes the 
distribution of the difference of the transverse momenta of the quarks px — 
(pi — P2)t/2 to be proportional to 1 + cos 2 9, where 9 is the angle between 
the vectors px and qx = (pi + P2)t/2. In the case of hadron interactions the 
polarizations of both gluons have to be taken into account as the correlation 
is more complicated. 

Although the situation considered here seems to be quite opposite to the 
parton model there is a certain limit in which our formulae can be trans- 
formed into parton model ones. Let us consider the pp case and assume 
now that the characteristic values of quarks' momenta p\x and p2x are much 
larger than the values of gluons' momenta qix, q?x 

< Pit >>< tor > , < Vrr >>< fer > (14) 

and one can keep only lowest powers of qix, qn- It means that we can put 
Qit — Q2T — everywhere in the matrix element M except the vertices. 
Introducing the polar coordinates 

d 2 QiT = ^dql T d9 1 (15) 

(and the same for q 2 x) and performing angular integration with the help of 
the formula 

p2tt 

/ d9iq^ T q\ T = nqfxS^ (16) 
J o 

we obtain 

[ 2n d9^ q ^- r d9 2 q -^ q -^M, u M Xa = 2vr 2 |l% \M part \\ (17) 
Jo y y Jo xx [xy) 

Here M part is just the matrix element in the parton model since the result 
is the same as that calculated for the real (mass shell) gluons and averaged 
over transverse polarizations. Then we obtain the cross section (||) in the 
form [1, 2]: 

da 

dy\dyld 2 p 1T 
|2 J_ f ^s{ql)f{y, ql T ) a a (g|)y(x, q\ T ) 2 2 
s) 2 J 4V2tt3 4V2tt 3 qiT q21 

a 2 (a 2 ) 

1 'iM^xGfaqfiyGfaqfi (18) 



2T 



(s) 2 

where s = xys is the mass square of QQ pair. 

On the other hand, the used assumption ( |T3D is not fulfilled in a more or 
less realistic case. The transverse momenta of produced quarks as well as 
gluon virtualities (QCD scale values) should be of the order of heavy quark 
masses. 
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3 RESULTS OF CALCULATIONS 



The functions <p(x, q 2 ) and <p(y, qf) are unknown at small values of q\ and q\ . 
Therefore we rewrite the integrals over q2 in eqs. (|3|) and in the form 

f°° 01 (do) 1 

/ dqj T 5(q 2 T ~ Pit ~ P2t) %i +2/2 - 1) 8 2 x ) \M\ 2 — = 

Jo ' qi V2 

\M\ 2 1 

4^2 7r 3 5(g 2T - pit - p 2T ) + y 2 - l)a,(Qg)a;G(ar, <2o)( J -2 L ) 92 -o— (19) 

<?2 2/2 

+ /, dql T 5(q 2T -p lT -p 2T )5( yi + y 2 -l)^f± V (qlx)\M\ 2 - 
J Ql Q2 V2 

where eq. (|l|) is used | and we use the same procedure for the integral over 
qi in eq. (5). 

In accordance with j|, || the value Ql = 4 GeV 2 is taken. The effective 
region of integration in the last integral of eq. (19) is restricted by the values 
of the matrix element which decreases fast at qf ^> uiq so the effective QCD 
scale is of the same order as in the parton model. 

Using eqs. (§), @, © an d eq. fll~9]) we perform the numerical calculation 
of heavy quark production. Two variants of gluon distribution are used. The 
first one is MT (S-DIS) distribution while the second is the result of ref. 
[|]. The last one was obtained using MT (S-DIS) boundary condition with 
inclusion of the shadow corrections changing the gluon distribution in the 
small x region. For heavy quark masses we used the values m c = 1.5 GeV 
and rrib = 4.7 GeV and the scale values are presented in eq.(19). 

In Fig. 3 we present the cross sections of charm and beauty production 
in pp collisions together with the lowest a s order parton model results with 
the same MT (S-DIS) gluon distribution. One can see that at moderate 
energies where not very small x values contribute, the difference in our results 
with standard MT (S-DIS) gluon distribution and with the distribution of 
ref. is negligible. An effect of screening corrections which decrease the 
growth of cross sections can be seen only in the case of charm production at 
LHC energies. Our cross sections are a little higher than the parton model 
predictions, especially in the case of beauty production and have more strong 
energy dependence. 

Some time at not too high energies, where the essential xi are not small 
enough the dashed curves (with the screening corrections) go higher than the 
solid ones. The explanation is that after the screening corrections were taken 

3 Please, don't mix the values of Q\ here and q${x) which was discussed in the 
Introduction. 
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into account in the evolution equation 0, the values of ip(x, q 2 ) were slightly 
changed in comparison with the original MT (S-DIS) set. 

Our predictions for the transverse momentum distributions of produced 
charm and beauty quarks are presented in Fig. 4. Here also we can see the 
screening effect for c-quark production at \fs =16 TeV. At large pt and not 
very high energies when comparatively large a;-values contribute there is a 
small difference in cross sections obtained with gluon distributions from || 
and ||. It is the consequence of a small difference in the values of (p(x, q 2 ). 

In fig. 5 the results for beauty production in pp collisions are compared 
with the parton model predictions and with the existing experimental data 
at \/s = 1800 GeV taken from ||. Two parton model curves are evaluated 
using MRSA [nj set of parton structure functions. They are shown as dash- 
dotted lines. The parameters m b = 4.5 GeV, ji = fi /2 (fi Q = ^Jm 2 + pj), 
A 5 = 300 MeV and m b = 5 GeV, /i = 2/i , A 5 = 151 MeV were used for the 
upper and lower dash-dotted curves, respectively. Both our curves (solid and 
dashed ones) calculated in LO are more smooth and pass between two curves 
of the parton model where NLO corrections was accounted for. Our curves 
are in good agreement with DO data and with CDF data starting from px > 
10 GeV/c. In all cases our approach results at high energies in slightly more 
broader p^-distributions m comparison with the parton model predictions 
using the same quark mass and scale values. 

Our predictions for rapidity (y) distributions of charm and beauty quarks 
produced in pp collisions at = 1800 GeV are presented in Fig. 6. Small 
screening effects can be seen at large values of \y\. 

The total cross sections of charm and beauty quark photoproduction are 
presented in Fig. 7 together with the parton model predictions in LO. Here 
the situation differs slightly from the case of pp collisions, our cross sections 
of charm production at low energies are smaller than the parton model pre- 
dictions. However the energy dependence of our cross sections is more strong. 
The screening corrections are negligibly small until the maximal energies of 
HERA. 

Our predictions for p T and rapidity dependences of charm and beauty 
quark photoproduction at y^i = 200 GeV are shown in Figs. 8 and 9. 

We calculated also the transversal and longitudinal cross sections of charm 
and beauty quark electroproduction in kinematical region of HERA. They are 
presented in terms of the functions F 2 (x, Q 2 ) in Fig. 10. The screening effects 
here are very small too. The ratios o"l/ot, presented in Table 1 increase with 
Q 2 (as for small Q 2 — > <jl should tend to zero) but never exceed the value 
0.4. 

To study the role of small q\ gluon virtuality we have changed the value of 
Q\ in eq. (19) from 4 GeV 2 to 1 GeV 2 . In this case (when the parton model 
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formulae were used only for q\ < 1 GeV 2 , while the true q 2 r dependence of the 
matrix element \M\ was taken into account beginning from q\ > 1 GeV 2 ) the 
DIS cross section (i.e. F 2 values) for charm electroproduction were increased 
by « 5 - 7% for W = 200 GeV and Q 2 = 4 - 64 GeV 2 . Thus the difference 
between our calculations and the parton model results increases. 

Let us note that even in the region of Q 2 = 16 - 64 GeV 2 that is much 
larger than the charm quark mass square m 2 Q = 2.25 GeV 2 we observe in Fig. 
10a strong scaling violation for the function F 2 . It reflect the fact that the 
essential transverse momenta of the produced quarks pxi are smaller than 
^fCp and up to Q 2 ~ 64 GeV 2 they stay comparable with ijiq value. Thus 
the scaling behaviour of F 2 will be restored only at very large values of Q 2 . 



4 CONCLUSION 

One can see that the accounting for the virtual nature of the interacting 
gluons as well as their transverse motion and different polarizations result 
in some qualitative differences with parton model predictions. However the 
numerical values of these differences at the considered energies are possibly 
smaller than the uncertainties of the predictions. Some screening effects 
are expected to be essential only in the case of charm production at LHC 
energies. 

Our calculations of heavy quarks production based on the lowest order 
matrix elements for the elementary subprocesses since our goal was to incor- 
porate the effects of gluon virtualities and polarizations which are important 
at small x. The results obtained show this effects to be actually existing. 
The high order corrections can be accounted for by using the effective value 
of -ftT-factor taken from the parton model calculations. As it is known the 
main (~ 0(iras)) part of the if-factor comes from the Sudakov double log 
formfactors, which are the same in our scheme as well as in the parton model. 
We can expect therefore our results are to be rescaled by approximately the 
same factor and the qualitative picture will be true also in this case. 

This work is supported by Russian Fund of Fundamental Research (95-2- 
03145). 
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Table 1 



The values of &l / °r ratios for the case of charm and beauty electroproduction 
at W = 200 GeV. 



QQ 


Q 2 = A GeV 2 


Q' 2 = 16 GeV 2 


Q2 = 64 GeV 2 


cc 


0.09 


0.21 


0.31 


bb 


0.012 


0.044 


0.13 
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Figure captions 



Fig. 1. Low order QCD diagrams for heavy quark photo- and electropro- 
duction via photon-gluon fusion. 

Fig. 2. Low order QCD diagrams for heavy quark production in pp 
collisions via gluon-gluon fusion. 

Fig. 3. Our predictions for heavy quark production cross sections in pp 
(pp) collisions together with LO order results of parton model. 

Fig. 4. The calculated px distributions of charm (a) and beauty (b) 
quarks produced in high energy pp (pp) collisions. 

Fig. 5. Our calculations of beauty production cross sections in pp (pp) 
collisions at ^/s = 1800 GeV (solid and dashed curves) together with NLO 
order results of parton model (dash-dotted curves) and experimental data of 
CDF (a) and DO (b) Collaborations. 

Fig. 6. Our predictions for rapidity distributions of heavy quark produced 
in pp (pp) collisions at y/s = 1800 GeV. 

Fig. 7. Our predictions for heavy quark photoproduction cross sections 
together with LO order results of parton model. 

Fig. 8. The calculated pr distributions of charm and beauty quarks 
produced in ^p interactions at y/s = 200 GeV. 

Fig. 9. Our results for rapidity distributions of heavy quark produced in 
7P interactions at \/s = 200 GeV. 

Fig. 10. Our predictions for functions F 2 (x, Q 2 ) in the cases of charm (a) 
and beauty (b) electroproduction in HERA kinematical region. 
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